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The possibility of selective harmonic excitation in a one-dimensional plasma that obeys the Vlasov 
equation is discussed. Electron beams, collimated in velocity space, are used for this purpose. Condi- 
tions on the velocity distribution function to achieve selective excitation are given. The velocity dis- 
tribution function is expanded in Gram-Charlier series and the expansion coefficients are given. Bounds 
on truncation errors of the expansion are derived and compared with computer results. 

Key words: Electron trapping; non-linear; numerical; selective harmonic excitation; two-stream; 
velocity distribution; Vlasov equation. 

Introduction 

This paper deals with the feasibility of a numerical study of higher harmonic excitation of the 
non-linear Vlasov plasma by two counterstreaming electron beams. In section 1 it is shown that 
an electron beam with a high degree of collimation in velocity space can couple to the phase 
velocity of one of the higher harmonics. The coupling should provide preferential excitation of 
the harmonic. The collimation is achieved by expressing the velocity distribution function in the 
form exp [— (v ±a) 2 /2e] with e being the collimation parameter that determines the halfwidth 
of the velocity distribution curve. The distribution function is expanded in Gram-Charlier series 
and the expansion coefficients are derived in section 2. The necessary transformations of Hermite 
polynomials are given in the paper "Temperature Renormalization in the Non-Linear Vlasov Prob- 
lem," immediately preceding this paper. In section 3 are discussed convergence properties of the 
Gram-Charlier expansion and a criterion is given that allows one to calculate an upper bound on 
the number of terms necessary to represent the initial distribution function to the desired degree 
of accuracy. 

Comparison of analytic results with those obtained on the computer are discussed in the 
appendix. 

1. Selective Excitation of Harmonics 

Numerical work on the solution of the non-linear Vlasov equation by the method of eigenfunc- 
tion expansion performed by one of the authors [1] l brought out some interesting facts associated 
with the non-linear character of the equation. One of these is that the ratio of plasma frequencies 
of successive harmonics is not an integer. For example for A; = 0.5 

a>2/a>i=1.86 (1.1) 



1 Figures in brackets indicate the literature references at the end of this paper. 
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where k is the wave number of the fundamental, k = 27r/L. The phase velocity of a standing wave 
in a plasma is given by 

(Op 

v» = -£ (1.2) 



with ojp being the plasma frequency. Since the ratio of A:'s in a periodic system is an integer, we can 
see from eqs (1.1) and (1.2) that the phase velocities of the fundamental and the second harmonic 
do not coincide. In the example quoted above the phase velocities of the fundamental and the 
second harmonic are respectively 

v P hi = 2.1l and %/,2 = 2.52. 

The difference in phase velocities is Af = 0.19. This fact suggests the possibility of a computer 
study of selective excitation of harmonics by two counterstreaming beams of electrons traveling 
through the plasma. The electron beams must be sufficiently well collimated in velocity space so 
as to couple to the desired harmonic. The amount of trapping as a function of the distance from the 
phase velocity was given by J. Dawson [4]. The functional dependence is rather complicated and 
will not be given here. It is simpler to use numerical results to examine the possibility of selective 
excitation of harmonics. Figure 1 shows f(v) that satisfies the non-linear Vlasov equation 2 for 
k = 0.5, a = 0.1 at the dimensionless time £ = 32. 

The amplitude of the second harmonic was at all times at least an order of magnitude less than 
that of the fundamental. The width of the shoulder as can be seen from the insert in figure 1 is 
about 1.54, with the phase velocities of the fundamental and the second harmonic located on this 
shoulder. Thus if an electron beam with a velocity distribution with a halfwidth of the order of 
separation of the two phase velocities were centered on the phase velocity of the second harmonic, 
the interaction of the second harmonic with the electron beam would excite this harmonic. The 
effect of the beam on the fundamental would be much less pronounced because of the lower density 
at the phase velocity of the fundamental. 

Since the initial distribution function [1] was chosen to be Maxwellian, a natural choice for the 
two-stream distribution function is also a Maxwellian plus a combination of Maxwellians centered 
at a and — a. 



f(v,t=0) = 



v2_ (v-a)2 (v+ap "I 

2 + ye 2€ + ye 2e (1.3) 



27r(l + 2yV^) 

where — 7= -j=- is the normalization constant. The parameter e determines the halfwidth 

V277-(l + 2yVe) 
of the two-stream distribution function. To obtain a halfwidth of At>, e must satisfy the relation 

(At;) 2 

^-1. (1-4) 

For Af = 0.19 this gives the value 

e=1.8lXlO- 2 . 



2 Unpublished results. 
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The individual components of the two-stream distribution functions are expressed in the form 

(»-«)2 J,2 



e *' =*2CnHe„(v)e> 



(v+a)2 



(1.5) 



e 2 « =^D n He n {v)e 2 



where //e„ are Hermite polynomials. The convergence of the expansion, the value of its coeffi- 
cients, and the number of terms in the expansion will be treated in detail below. 
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FIGURE 1. f(v) versus v for k=0.5 at t = 32. 

The dashed line gives e _l ' 2/2 /V27r for comparison. The insert represents/(w)/(e- ,,2/2 /\/27r). 
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2. Expansion Coefficients for the Two-Stream Problem 

We want to express one of the side peaks of the two-stream velocity distribution function of 
eq (1.5) in the following form: 

(r-q)2 _ y£ 

e 2 < =g(v)e 2 (2.1) 

where 

g(i>) = J CnHe n (v). (2.2) 

rc = 

Coefficients for the other side peak can be obtained immediately by substituting a for — a, as can 
be easily seen by comparing eqs TR (1.14) 3 and TR (1.15). The expansion coefficients, C n , are 
given by the integral: 

1 f« -- 

C n = ,— g(v)e 2 He n (v)dv. (2.3) 

\J2irn\ J -oo 

Substituting the expression for g(v) from eq'(2.1), we obtain 

1 foe (v-a) 2 

C n = — =— e 2 * He n (v)dv. (2.4) 

V2rrnl J-oo 

Let us now make the transformation 
This leads to the integral 

Here fi = e~ l was introduced to express the argument oi He n in the form given in eq TR (1.15), 

and a = — f= • Since He = \ the integrad in eq (2.6) can be considered as a product 

Ve 

He n = He*He n . (2.7) 

//e n f — -pr j is expressed in terms of ^He n (v) by eq TR (1.15). From the orthogonality conditions 
we see that the only term in the sum that contributes to the integral is the term with p = n. 
Substituting this term for He n ( — j= I, integration can be carried out and we obtain: 

C n =i €l / 2 (l-e)»/^ e „(^=)- (2.8) 



3 TR-refers to "Temperature Renormalization in the Non-Linear Vlasov Problem" by the same authors, immediately preceding this paper. 
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Substituting the values of C n obtained in eq (2.8) into eq (2.2), we can now express eq (1.3) in terms 
of its Hermite polynomial expansion. 



f(v,t = 0)=- 



r-l-l 



/ 2^(l+2yV7) L 



1 + y Z, i He " I 77== - ) He n (v) 



VT^ 



+y g*^ H<i ^ )H , M] . 



(2.9) 



Since Hermite polynomials with odd indices are odd functions, eq (2.9) reduces to its final form 



-vm 



€ l/2(l- € )« 



V27r(l + 2yV€) L So ( 2n ) ! Wl-e 

3. Convergence Properties of the Expansion 



He 2n (v) 



(2.10) 



In numerical work it is desirable to truncate the expansion after the smallest number of terms 
for a given truncation error. Using the asymptotic form for Hermite polynomials [6] (see TR appendix 
I for the discussion) 



„ , . e^ 4 r(nH-i) r ( m/2 x i \, n/ 1/2 . 

where N=2n-\-l. For large n, we obtain the asymptotic expression for 

e l u- € )i s?4-*L r 2 (^i + D 

r //p (^ ~ — — P 4 4(l-e) — 

L„tie H (v) n[ 2 - n e T 2 (^ + l) 



(3.1) 



(3.2) 



with C„ defined by eq (2.6). 

In the above expression the term 0(n~ 1/2 ) has been neglected and cos ( N l/2 — ^= — - rnrj has 

been replaced by unity. 

Using the well-known relation between the r -function and the factorial 

r(z + l)=z! 

we obtain the asymptotic expansion for the terms in eq (2.10) 

e 1 ' 2 (2n) ! 



(3.3) 



C 2n He 2n {v) ~ € K 



2 2 » (nl) 2 



(1-6)" 



(3.4) 



yhere K = — + 



V L , (T 



4 4(1-6) 



. Using Stirling's formula for the factorials, we finally get: 



C 2n He 2n 



e m e K (i_ e )« 



tt v n 



(3.5) 



The asymptotic expression for C>„He>„ in eq (3.5) tends monotonically to zero with increasing n. 
This can be shown easily by taking the logarithm of C>„Hem 



In CzitHezti ~ In B + n In (1— e) — £ In n 
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(3.6) 



362-120 O - 69 - 2 



wher 



/2 P K 



e 1/z e 



Vtt 



Since < e ^ 1, In (1 — e) is a negative number and the logarithmic terms in eq (3.6) go to — °°. 
Thus lim C 2n He 2n = 0. (3.7) 

Let us now find the truncation error introduced into the function by neglecting all the terms 
in eq (2.10) past the term He 2 N- This error can be represented in the form: 



yhere 



"-■ •= * •"('»+2)T' """• ( vfa) H -" w 1° -"' 



<9 = 2/V2^ (l + 2yV€). 



(3.8) 



Expressing all He 2 N+ 2 +p, p = to ™ in their asymptotic form we see that since 1 — e < 1, the geo- 
metric progression in eq (3.8) can be summed and we finally obtain a bound on the truncation error 



i+- 



Rn+i < 



,4 4(l-€) 



(l-e) J 



Vtt e" 2 (/V+l)" 2 



(3.9) 



"i i rn i i i i i 1 1 



I I I LI 




10 



FIGURE 2. N versus e for two-stream function. 
Computer results and analytic results for = 3.0, a = 2.0 
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The inequality in eq (3.9) makes it possible to calculate N for a given truncation error. This is 
useful in representing the initial two-stream distribution function in the Vlasov equation. Figure 2 
gives the number of terms in the expansion of the expression given in eq (2.10) as a function of e 
for a = 2.0 and v = 3.0 for a truncation error of R n < 10~ 5 ,S (see appendix). A graph of the number 
of terms obtained on the computer is included for comparison. The values of N for e = 0.9 and 0.8 
do not lie on the curve, which is not surprising since the asymptotic expression for He n is not 
very accurate for /V ~ 10. As A^ increases the points form a curve which bounds from above the 
computed curve. This again is not surprising, since the approximations made in eqs (3.2) and (3.8) 
made the truncation error somewhat bigger than it actually is. 



4. Appendix 

A computer program was written in Fortran to calculate the terms of the expansion oif(v, t = 0) 
of eq (2.10). The He 2n (v) and C 2n were calculated recursively. Scaled polynomials, He 2 n, were used 
to stay within the range of the computer in single precision. 



He 2n (v)=- 7 =He 2n (v) (A.l) 

V (2n) ! 



C 2 »=V(2n)! C 2 „. 



A partial sum, P n , of the first n terms was formed and more terms were added until the relative 
error ratio test of eq (A.2) was satisfied for one hundred successive terms. 

|R.E.(P H )| = |P y~ Pa l <10-». (A.2) 



Using eqs (2.10), (3.5), (3.9), and (A.2), the relation of this ratio test to the truncation error is given 
by eq. (A.3). 

6 e A (l-e) v+l 
R N+1 < ^ , } -<10-*S (A.3) 

e" 2 V^ (/V+l)" 2 
where 

OJe" 212 
e 2 

and/ is the un-normalized function calculated on the computer. Figure 2 contains the computer 
results for t; = 3.0, a = 2.0 and a comparison analytic curve based on eq (A.3). Other values of the 
parameters investigated on the computer were 

6= 1.00, 0.30, 0.10, 0.07, 0.03, 0.01 
a =1.00, 2.00, 3.00, 4.00, 5.00 
v= 0.00 to 6.00 in intervals of 0.25 
y=0.1. 

Figure 3 contains a set of results for t> = 3.0, a =1.0, 3.0, 5.0. 
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Figure 3. N versus efor two-stream function. 

Computer results for v — 3.0 
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